Entanglement renormalization of anisotropic XY model 
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The renormalization group flows of the one-dimensional anisotropic XY model and quantum Ising 
model under a transverse field are obtained by different multiscale entanglement renormalization 
ansatz schemes. It is shown that the optimized disentangler removes the short-range entanglement 
by rotating the system in the parameter space spanned by the anisotropy and the magnetic field. 
It is understood from the study that the disentangler reduces the entanglement by mapping the 
system to another one in the same universality class but with smaller short range entanglement. 
The phase boundary and corresponding critical exponents are calculated using different schemes 
with different block sizes, look-ahead steps and truncation dimensions. It is shown that larger 
truncation dimension leads to more accurate results and that using larger block size or look-ahead 
step improve the overall calculation consistency. 

PACS numbers: 05.10.Cc, 75.10.Pq, 02.70.-c, 03.67.-a 



The real-space renormalization-group (RSRG)fi re- 
volving around the coarse-graining and rescaling trans- 
formation, has been proven to be a very useful tool in 
the understanding of the critical phenomenons and the 
quantum many-body system, whose difficulty lies in the 
exponential growing of the Hilbert space with the system 
size. The original RSRG, introduced by Wilson based 
on the block spin idea^ addresses the problem by di- 
viding the system into blocks and truncating the block 
Hilbert space to a subspace spanned by a few eigenstates 
of the lowest energies i Later, White suggested that one 
should consider the interplay between the block and its 
environment in the truncating and promoted the den- 
sity matrix renormalization group (DMRG)^ In DMRG 
algorithm, the states to be retained are the eigenstates 
with the largest eigenvalues of the reduced density matrix 
of the ground state of the block and its environment. 3 
Since the reduced density matrix measures the entan- 
glement between the block and its environment, it is 
later understood that the performance of DMRG de- 
pends on the entanglement in the ground stated Un- 
fortunately, near the quantum critical point, the ground 
state has large entanglement and hence one needs large 
truncated dimension to get accurate result. To solve 
this problem, Vidal proposed a new entanglement renor- 
malization method, multiscale entanglement renormal- 
ization ansatz (MERA)j^S by introducing an additional 
unitary transformation U, the disentangler, which acts 
on the boundary of the adjacent blocks to remove the 
short-range entanglement (SRE) before performing the 
coarse-graining. MERA has been shown to be a suc- 
cessful numerical scheme in a lot of different physics 
systems, such as one-dimensional^ and two-dimensional 
quantum spin system^ interacting Fermions^ boundary 
critical phenomenal However, there are still some very 
important problems about MERA that need to be an- 
swered, such as how exactly does the disentangler re- 
move SRE, what is the difference among the different 
MERA scheme? In this paper, we try to understand these 
problem by applying MERA to study the anisotropic XY 



model (AXY) under a transverse field. 

The Hamiltonian of AXY is defined as H 
within 



E^. 



M+li 



Hi, i+ x = ~[(l+l)S?Sf +l +(l-T)S?S? +1 ]-~(Si+S? +1 ), 

where i is the site index, S a = a a /2 (a = x, y, z) are the 
spin components represented by the Pauli matrix a a . J, 
7 and h stand for the interaction strength, anisotropy of 
the interaction and transverse field respectively. In the 
limits of 7 = and 1, the model becomes the XY model 
and the Ising model in a transverse field (ITF) respec- 
tively. The model is exact soluble by using the Jordan- 
Wigner and Fourier transformations. 11 The ground state 
of AXY in the regime < 7 < 1 belongs to the quan- 
tum Ising model universality class yti The system exhibits 
three critical lines at x c (-f) — h c ("f)/ J = ±1 and 7 C = 0. 
Under weak magnetic field (|a;| = \h/J\ < 1), the system 
is in the ferromagnetic phase. When the field increases 
to the critical value x — x c (j) = 1, the system under- 
goes a quantum phase transition (QPT) then turns into 
the paramagnetic phase when the field further increases 
above the critical value. Due to the exact solubility and 
rich physics, AXY and its special case ITF have been 
extensively studied to understand the nature of QPT, 
especially the role of the entanglement in QPTJ^r— It 
also provides a test field for new numerical schemesi 15 i 16 
The phase diagrams and the critical exponents of AXY 
and ITF have been obtained using various RSRG and 
DMRG schemes. In this paper, we study these models 
by using different MERA schemes. 

MERA can be understood as a quantum circuit or 
renormalization group (RG) transformation! 5 i 6 i 17 Here 
we focus on the RG point of view. As demonstrated 
in Fig. [TJ the system is divided into blocks of n sites, 
the RG transformations are performed by applying a se- 
rial of discntanglers U and isometries W on these blocks. 
The coarse-graining is implemented by the isometry W, 
which maps the Hilbert space on a n-site block into a new 
space on a coarse-grained site. The new Hilbert space is 



truncated to a subspace of dimension \ small enough 
for one to carry out the calculation, but large enough to 
represent the system faithfully. The purpose of the dis- 
cntangler is to remove SRE so that one can use small \ 
to get accurate result even in the critical regime. The 
applying of the disentangler U and the isometry W lifts 
the system on the original sites to a new system on the 
coarse-grained sites. By comparing the properties of the 
Hamiltonian or the ground state of the original system 
and the corresponding coarse-grained version, one can 
get the renormalized parameters that define the coarse- 
grained system. Repeating the process one then gets a 
well defined RG flow. The RG flow of the entanglement 
was obtained for the ground state of QIT and AXY mod- 
els using a special MERA scheme targeting at the max 
entangled states by removing the unentangled modes. 5 
It was shown that the disentangler indeed reduces the 
inter-block entanglement. However, this special MERA 
scheme is designed to obtain max entangled states but 
not suitable for the general purpose. Moreover, the RG 
flow in the traditional sense, i.e. the change of the physics 
parameters that define the system Hamiltonian under the 
RG transformation, and how the disentangler affects the 
flow remain to be further studied. 
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FIG. 1. MERA scheme 



In order to target the ground state, we assume that 
the ground state takes the form of the block mean field 

state 18 after r-layer coarse-graining. 
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where \<jr- }, whose form does not depend on j, is a wave- 
function defined on the j'-th block at the r-th coarse- 
grained layer. In MERA's language, we use translational 
invariant MERA, but set the truncation dimension of 
the top (r + 1-th) layer isometry to be 1. The disen- 
tanglers and isometries corresponding to this block mean 
field state are optimized by the algorithm for the trans- 
lational invariant MERA presented in Ref. |fl. For the 
coarse-grained Hamiltonian to have same symmetries as 
the original one, the disentanglers and the isometries 
are chosen to preserve the spatial reflection symmetry 
and global parity symmetry^ It should be noted that 
even though we might use more than one layer of coarse- 
graining to obtain the ground state in the form of Eq. (|2|) , 
the RG flow is solely determined by the disentangler U 
and W at the first layer. Our RG scheme is similar to 
the Born-Oppenheimer RG scheme used in Ref. [l8| with 



multiple look-ahead steps. However, we do not need 
to introduce some artificial parameters to describe the 
slow mode, which is nothing but the effect of the inter- 
block interaction and is naturally and self-consistently 
taken into account in MERA scheme. Like the origi- 
nal Born-Oppenheimer RG scheme, r can be viewed as 
look-ahead step. The scale invariant MERA-£ used to 
study the critical system can be seen as an approxima- 
tion of infinite look-ahead steps. In the following we will 
discuss the results for different MERA schemes, such as 
with/without disentangler, different block size n, look- 
ahead step t. For simplicity, the schemes are denoted by 
a string 'nrxZ', where n is the block size, r is the look- 
ahead step, x can be 1 for the simple block mean field 
state or oo for the scale invariant scheme. Z is 'D' for 
schemes with or T for without disentangler respectively. 
We first focus on the simplest case when only two 
states, one with even parity and the other with odd par- 
ity, are retained so that we can track the RG flow ana- 
lytically. Under the restriction of symmetry requirement, 
the disentangler takes the form of U = Ue + Uq with Ue 
and Uo apply on the even and odd parity subspace re- 
spectively. For the disentangler acts on sites and 1, Ue 
and Uq take the following two possible forms 



os 8 f-i 



COS V ( z 






Uo = (1 - (tfof )/2 or ffiof + a y af)/2 . 



(3) 

(4) 



The disentangler keeps the form of the AXY Hamiltonian 
Hii + \ on the boundary sites unchanged, but maps the 
anisotropy 7 and magnetic field x = h/J to 7' and x' , 
with 



7' = 7 cos 29 + x sin 28, 
x = —7 sin 26 + x cos 29. 



(5) 
(6) 



That is, it rotates the vector (7, x) in the parameter space 
by 29. Since SRE depends on these parameters; 12 ' 20 it 
is possible for the disentangler to remove the inter-block 
entanglement by suitable rotation. It is noted that even 
if one starts from ITF model in which 7 = 1, applying 
of disentangler renders it to general AXY model with 
7^1 unless one is limited to 9 — 0, which is equivalent 
to without applying disentangler. This is quite different 
from the previous RSRG schemes for ITF where the RG 
transformations do not change the form of ITF Hamilto- 
nian. The change from ITF to AXY should not affect the 
calculation of the critical properties since the anisotropy 
is irrelevant term and AXY belongs to same universality 
class for for all regime < 7 < f. From this aspect, 
MERA actually provides a possible way for one to iden- 
tify the universality class. 

To see the effect of the disentangler on SRE, we use 
the concurrence of the ground state of the two qubit 
system 20 composed of the spins on the two adjacent sites. 
It is known that the concurrence of the two qubit sys- 
tem stays at f for r\ = ^/■y 2 + x 2 < 1 and decreases as 
7/77 when r\ is larger than one. At exactly r) = 1, the 




FIG. 2. (Color online) Concurrence of ground state of the 
two qubit AXY system at different RG stages as functions 
of magnetic field under different anisotropies: Dashed (red) 
curves: Original system; Solid (black) curves: After applying 
the optimal disentangler; Dotted (blue) curves: After one RG 
transformation. 



concurrence is 1 — 7/2 which is a local minimal when 
x < 7 3 (4 — 7)/(2 — 7) 2 . 20 Strictly speaking, this concur- 
rence can only measure the entanglement of the isolated 
two qubit system and should not be regarded as the ex- 
act SRE between the two boundary spins. Nevertheless, 
it provides some qualitative information of SRE and the 
simple form enables one to understand the role of disen- 
tangler analytically. In Fig. [21 we plot the concurrence of 
the ground state of the two qubit system after the apply- 
ing of the disentangler for different magnetic fields and 
anisotropies. For comparison, we also plot the concur- 
rence of the original and coarse-grained two-qubit sys- 
tem. The figure clear shows that the disentangler indeed 
reduces SRE. Analytically, the optimized 9 is larger than 
zero when 77 < 1 and continually approaches to negative 
values when 77 becomes larger than 1. When 77 = 1, the 
optimized 6 is about 0. One can see from these results 
that, the optimized disentangler reduces SRE by map- 
ping the system to a new one of the same universality 
class but with smaller entanglement. 




FIG. 3. (Color online) RG flows and phase boundaries 
obtained from different MERA schemes: Solid black lines 
with arrow, RG flows from 321D scheme; Solid red curve, 
phase boundary with 321D scheme; Dashed blue curve, phase 
boundary with 52 ID scheme. The critical fixed points are 
marked by cross. 

In Fig. [31 we present the flows in the parameter space 



under RG transformation and the phase boundary de- 
fined by the flows. It is seen that the system has three 
critical lines at 7 = 0,.t < 1, and < 7 < l,x = ScCt) 
that divide the parameter space — 1 < 7 < 1 and 
— 00 < x < 00 into four parts (since the diagram is sym- 
metrical under x — > —x and 7 — > —7, only one part is 
shown in the figure). Starting from 7 > 0,x < x e (7), 
RG transformations eventually bring the system to the 
attractive point 7 = l,x = which represents the fer- 
romagnetic phase. While staring from 7 > 0, x > x c (j), 
it is brought to another fixed point at 7 = 0, x = 00, 
corresponding to the paramagnetic phase. On the crit- 
ical line x = x c (j), there is another critical fixed point 
< 7 = 7 C < 1. In Table. [J we list the critical field for 
7 ~ and 1, critical fixed point and some of the criti- 
cal exponents obtained from 321D scheme. To show the 
effect of the disentangler, we also list the corresponding 
results from 3211 scheme, which is identical to 321D ex- 
cept that it does not have disentangler, for comparison. 
One can easy see that adding the disentangler greatly 
improves the accuracy of the phase boundary and the 
critical exponents. For example, without the disentan- 
gler, the critical field for Ising model is x c (l) = 1.291 
from 3211 scheme while with disentangler it is 1.081. 
Moreover, these two schemes give different critical fixed 
points. For 3211 schemes, it is j c = l,x c = 1.291 which 
is in consistence with the results from standard RSRG 
and DMRG schemes. However, with disentangler the 
critical fixed point becomes j c — 0.801, x c — 1.071 for 
32 ID scheme, clearly deviates from that of 3211 scheme. 
The difference in the positions of critical fixed point is 
a result purely from the disentangler. It also reveals the 
role of the disentangler in MERA scheme. Namely, the 
disentangler changes the position the critical fixed point 
so that the critical system has smaller SRE and can be 
better represented by small truncation dimension. 

The number of the look-ahead step also affects the ac- 
curacy. Generally speaking, larger r leads to more accu- 
rate results but also larger calculation cost. The optimal 
choice of r can be understood from RG and MERA point 
of view. The wave-function obtained here is a variational 
one in the product state form. For system that is not 
close to the critical regime, a few RG transformations 
should bring the system to the attractive fixed points 
such as the paramagnetic or ferromagnetic states, whose 
ground states are the simple product states. Therefore 
for system away from critical regime, one can use the 
wave-function in the form of Eq. @ with small r to rep- 
resent the ground state accurately. Further increasing 
r should only have small effect since the disentanglers 
and isometries on the top layers would be identical to 
each other as the system flows to the trivial fixed points. 
As the system approaching the critical regime, more RG 
transformations, hence larger t, are need to bring the sys- 
tem to the product state. Exactly at the critical points, 
one needs infinite r to represent the ground state faith- 
fully using the product state form. In this case, the disen- 
tanglers and isometries on the lowest layers are different 



from each other since the irrelevant terms are different at 
different RG stages. As one goes up to the higher layer, 
these irrelevant terms eventually vanish and the system 
is brought to the critical fixed point. After that, the dis- 
cntanglers and isometries will be the same for different 
layers due to the scale invariance. In practice, one can 
choose a few "free" layers with different disentanglers and 
isometries then put scale invariant layers on top of the 
"free" layers. This scale invariant MERA scheme was 
used to study the properties of various critical systems^ 
It was believed that the scale invariant MERA can only 
be applied to the critical system. However, we argue that 
since it can be seen as an approximation of MERA with 
infinite look-ahead steps, one can expect that the scale 
invariant MERA can represent the system faithfully in 
the vicinity of the fixed points not matter they are criti- 
cal or noncritical. 

The results with different look-ahead steps are listed in 
TableQ] At first glance, it seems that different look-ahead 
steps do not lead to any significant differences in phase 
boundary and critical exponents. Product state MERA 
and scale invariant MERA also only give slightly differ- 
ence in the phase diagram and critical exponents. How- 
ever, further studies show that larger look-ahead steps 
and scale invariant scheme improve the accuracy on the 
other physics quantities, such as the susceptibility, and 
the consistency of these quantities. For example, the 
susceptibility, proportional to —d 2 E g (h)/dh 2 , should di- 
verge at the phase transition point. However, with 321D 
scheme, the susceptibility of ITF merely has a smooth 
peak at x = 0.95, which is far away from the critical 
field determined by the RG flow. As one increases the 
look-ahead step, the peak sharpens and moves towards 
the critical field. Using scale invariant MERA further 
improves the consistency. Using 33ooD scheme, the po- 
sition of the peak already consists with the critical field. 
Further increasing of "free" layer number further sharp- 
ens the peak but has little effect on the position. 

We now turn to the accuracy of RG for different block 
size. In Table. Q] we list the critical fields and critical ex- 
ponents obtained from different block sizes. It is known 
that the accuracy of the traditional RSRG can be im- 
proved by increasing the block sizei 15 ' 16 i 18 ' 21 ' 22 This fact 
is justified by the results without disentangler. As shown 
in the table, x c (l) is improved from 1.291 of the 3-site 
block scheme to 1.233 of the 5-site block scheme. How- 
ever with the disentangler, increasing the block size leads 
to slightly worse results, at least when the block size is 
small. Both the phase boundary and the critical expo- 
nents are getting a bit worse off when the block site in- 
creases from 3 to 5 and further to 7. However, like in- 
creasing the look-ahead step, increasing the block size 
also improves the consistency of the overall calculations. 
The 51ooD and 52ooD schemes product the susceptibility 
peaks similar to that from 32ooD and 33ooD schemes re- 
spectively. From this aspect, the increasing of block size 
does improve the overall accuracy in some sense. Follow- 
ing the logics of Ref. [22J, without the disentangler, the 



accuracies of the critical field and the critical exponents 
increase asymptotically as inverse of the block size or log- 
arithm of the block size. Since the results obtained from 
MERA with disentangler are always better than that 
without, it is expected that the accuracy on the phase 
boundary and critical exponents can be improved by in- 
creasing the block size when the block is large enough. 
However, increasing block size to improve the accuracy of 
the phase boundary and critical exponents is not a good 
strategy due to the following two reasons. One is the 
non-monotonic dependence of the accuracy on the block 
size and the formidable calculation cost of MERA scheme 
with larger block size. The other reason is that, as the 
block size increases, the disentangler becomes close to 
identity and the difference between the results with and 
without disentangler gradually diminishes. This trend 
can be seen by comparing the results from 321D, 3211, 
521D and 5211 schemes listed in Table, fl] In practice, use 
small block size but larger truncation dimension is a more 
feasible approach for one dimensional system with short 
range interactions. The optimal scheme for this kind of 
system is ternary MERA, but the other block sizes are 
also acceptable when desired. 
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0.125 


0.194 0.261 
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TABLE I. Phase boundary and critical exponents under dif- 
ferent MERA schemes. The scheme is denoted by nrxZ, 
where n is the block size, r is the look-ahead step, x can be 1 
for the simple product state or oo for scale invariant scheme, 
and Z is 'D' for schemes with or T for without disentangler 
respectively. The additional '(4)' in the last column means 
that the results are obtained with 4 states retained. 

We now turn to RG with more than 2 states retained on 
each coarse-grained sites. Theoretically, when one keeps 
2 l states, one can treat the coarse-grained site as a multi- 
site composed of I sites and write down the corresponding 
renormalized Hamiltonianj 16 i 23 However, the MERA RG 
transformation make the renormalized Hamiltonian very 
complicated even with only 4 states retained. Besides the 
next nearest neighbor two-body interaction resulted from 
the normal RSRG process, MERA also introduces addi- 
tional irrelevant two-body terms (S 1 ? — Sf +1 ) and SfS? +1 , 
three-body terms like Sf Sf +1 Sf +2 and four-body terms 
such as S*S? +1 Sf +2 S* +3 . Tracking the RG flow for AXY 
is a tedious task even for the 4-state case. A simpler 
way to draw the phase boundary is to look at the sus- 
ceptibility, whose peak position consists with the phase 
boundary when the look-ahead step is large enough. The 
critical field for 7 = 1 and the corresponding critical ex- 



ponents are calculated and listed in Table. U for \ = 4. 
One can see that using larger truncate dimension indeed 
greatly increases the accuracy of phase boundary and the 
critical exponents as it should have been. One can ob- 
tain very accurate critical exponents by using even larger 
truncation dimension.— 



study how the block size, look-ahead step and truncation 
dimension affect the accuracy. It is shown that increas- 
ing the block size and look-ahead step improve the over- 
all calculation consistency. Larger truncation dimension 
leads to much more accurate results in the phase diagram 
and critical exponents. 



In conclusion, we study the RG flow of one-dimensional 
AXY model and obtain the phase boundary and the 
corresponding critical exponents using different MERA 
schemes with the different block size, look-ahead step, 
and truncation dimension. It is understood that the dis- 
cntangler reduces the short range entanglement by chang- 
ing the system to a new system of the same universality 
class but with smaller short range entanglement. Espe- 
cially for AXY model, it is shown analytically that the 
optimized disentangler reduces the entanglement by ro- 
tating the system in the parameter space spanned by 
the anisotropy and the magnetic field field. We further 
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